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Abstract 

Due to the effects of surface electromigration, waves can propagate over the free surface of a 
current-carrying metalhc or semiconducting film of thickness /iq- In this paper, waves of finite am- 
phtude, and slow modulations of these waves, are studied. Periodic wave trains of finite amplitude 
are found, as well as their dispersion relation. If the film material is isotropic, a wave train with 
wavelength A is unstable if X/Hq < 3.9027. . ., and is otherwise marginally stable. The equation of 
motion for slow modulations of a finite amplitude, periodic wave train is shown to be the nonlinear 
Schrodinger equation. As a result, envelope solitons can travel over the film's surface. 
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I. Introduction 

Waves propagating over the free surface of a layer of water subject to gravity have been 
studied extensively for well over a century. In his ground breaking paper of 1847, Stokes 
showed that steady periodic wave trains of finite amplitude are possible, and that their 
dispersion relation depends upon their amplitude In these wave trains, the effects of 
dispersion and nonlinearity balance, yielding a steady wave form. 

In the 1960's, it was discovered that if the water is of finite depth ho and the wavelength 
A = 271 /k is sufficiently short, the wave train is unstable. Benjamin and Feir showed 
that the Stokes wave train is unstable to small modulational or side band perturbations if 
kho > 1.363 . . ., and this been confirmed experimentally 0, ^, ||, H. 

Benjamin and Feir's analysis holds only at early times, just after the instability has 
begun to develop. Zakharov and later Hasimoto and Ono @] derived an equation of 
motion for the envelope of a packet of finite-amplitude gravity waves which remains valid 
long after the onset of the instability. Their equation — the nonlinear Schrodinger equation 
— admits envelope soliton solutions that are intimately related to the ultimate state of the 
fiuid surface. 

Although it is not yet widely appreciated, there are nontrivial electrical free boundary 
problems. When an electrical current passes through a piece of solid metal or semiconductor, 
collisions between the conduction electrons and the atoms at the surface lead to drift of these 
atoms. This phenomenon, which is known as surface electromigration (SEM), can cause the 
surface of a solid to move and deform |, 0, 0, ||, |l|, 0, ||, |T|, P, |l|, [T|, H, |2l|]. The 



free surface of a conducting film moves in response to the electrical current fiowing through 
the bulk of the film, in much the same way that fiow in the bulk of a fiuid affects the 
motion of its surface. However, the analogy is not perfect — the boundary conditions are 
very different in the two problems. In addition, metals and semiconductors are anisotropic 
materials, and this distinguishes them from fluids as well. 

Waves can propagate over the free surface of a current- carrying fllm deposited on an 
insulating substrate, at least in the limit of small amplitude. These waves, whose propagation 
is induced by surface electromigration, are somewhat analogous to gravity waves, and will 
be called "SEM waves." The linear dispersion relation for these waves has been known for 
some time |TB[, and will be discussed briefly in Sec. III. 

It is natural to ask to what extent the nonlinear behavior of gravity waves is paralleled by 
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that of SEM waves. To address this question, in this paper I will study SEM waves of finite 
amplitude, and slow modulations of these waves. Periodic wave trains of finite amplitude 
(the analog of the Stokes wave) are found, as well as their dispersion relation. If the film 
material is isotropic, the wave trains are unstable for khQ > 1.609927. . . and are otherwise 
marginally stable. The equation of motion for slow modulations of the finite amplitude, 
periodic wave train is shown to be the nonlinear Schrodinger equation. 

It is worth mentioning that SEM is not just of academic interest: In modern integrated 
circuits, the current-carrying metal lines or "interconnects" carry very high current densities, 
and electrical failure can occur due to the effects of SEM. In particular, SEM can cause a 
small perturbation at the edge of a metal strip to become a slit-shaped void which propagates 
across the line until its tip contacts the opposite side of the strip; electrical failure is then 



complete 0, gl|,|2|, |23|. 
II. Equations of Motion 

Consider a single-crystal metallic or semiconducting film of uniform thickness ho deposited 
on the plane surface of an insulating substrate. We take the z axis to be normal to the 
film-vacuum interface and locate the origin O in this plane. (We will occasionally find it 
convenient to use another set of Cartesian coordinates {x', z') with x' = z' = z + and 
origin O' at the substrate surface.) The film's surface z = Q will be assumed to be a low 
index crystal plane. A constant electrical current flows through the film in the x direction, 
and the electric field within the film is Eq = Eqx. 

Now suppose that the upper surface of the film is perturbed. Let the outward-pointing 
unit normal to this surface be n. For simplicity, we shall restrict our attention to perturba- 
tions whose form does not depend upon y, so that the height of the film's surface above the 
substrate h = Hq + ( depends only on x and t (Fig. p. The upper film surface will evolve 
in the course of time due to the effects of SEM and surface self-diffusion. We assume that 
the current flowing through the film is held fixed. 

Clearly, the problem is two-dimensional (2D), and the dependence of all quantities on y 
will therefore be suppressed. The electrical potential $ = ^{x,z,t) satisfies the 2D Laplace 
equation 

V^<I) = (1) 
and is subject to the boundary condition n ■ V$ = on the upper surface and z ■ V$ = 
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FIG. 1: The current-carrying metal thin film. The height of the free surface above the substrate, 
/i, depends only on x and t. The outward-pointing unit normal to the free surface is n, and the 
electric field far from the perturbation is £"0- 



on the lower. More explicitly, we have 



(2) 



and 



$,(x,-/io,t) = 0, (3) 
where fx = df/dx and so forth. If the initial perturbation is localized, we also have 



(4) 



for X ^ ±00 and, furthermore, 

V^x,z,t)^-Eo (5) 

for X — i> ±cx) and —Hq < z < 0. 

We assume that the atomic mobility is negligible at the film-insulator interface, so that 
the form of that interface remains planar for all time. Further, in the interest of simplicity, 
we assume that the applied current is high enough that the effects of SEM are much more 
important than those of capillarity. The surface atomic current Ja is then proportional 



to the electrical current at the surface. Explicitly, Ja = — z/MgV$, where the v is the 
areal density of the mobile surface atoms, q is their effective charge, and M, the adatom 
mobility, in general depends on the surface slope C,x- When there is a net influx of atoms 
into a small surface element, it will move. The normal velocity of the film's free surface is 
Vn = qi'^ldsiMds^), where Q is the atomic volume and s is the arc length along the surface. 
Therefore 



Ct = qvVtd^ 



9,$(x,C,t) 



(6) 



where is the total derivative with respect to x. (Since C depends on x, we have 
= ^x^xXii) + Cx^z{.xX,i)i for example.) Together, Eqs. - (||) completely 
describe the nonlinear dynamics of the film surface. 
III. Finite Amplitude Periodic Wave Trains 

We wish to study the propagation of a disturbance whose amplitude is small. To this 
end, we put $ = —Eqx + and presume that ( and are of order e. 

We shall suppose that the crystal structure is invariant under the reflection x — ^ —x, so 
that M is an even function of the slope Cx- As a result, M{(x) = + c(^) + 0{e^), where 
Mq = M(0) and c is a dimensionless constant which vanishes for an isotropic material. 

If we work to order e, the equations of motion are linearized and we find that sinusoidal 
waves of the form ( = 2acos{kx — ujQt) propagate over the surface, c^o = ^o(^) is given by 
the linear dispersion relation 

= -PE,-, (7) 
a 



where /? = qM^uVL and a = ta.nh.kho The corresponding group velocity is 

(8) 



Vg = VO 



where Vq = —/3Eo/hQ is the phase velocity of waves of infinitely long wavelength. If the am- 
plitude 2a is not a constant but instead varies slowly with position, the resulting amplitude 
modulation will propagate with the group velocity. 

The goal of this paper is to go beyond the linear approximation and study wave trains 
of finite amplitude. To begin, we will find the analog of the Stokes wave propagating over 
deep water: we will consider the limit oo and restrict our attention to periodic wave 

trains with constant amplitude a. These simplifying assumptions will be relaxed in the next 
section. 
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Expanding ( and in powers of the amplitude, we have 

C=[a£ + H2a^£'^ + fi^a^S^ + ...]+ c.c. (9) 

and 

= [uiae^'S + V2a^e^^'£^ + v^a^'^^'E^ + v^a^e^'E + ...] + c.c, (10) 

where 8 = e*^'^^"'^*), "c.c." denotes the complex conjugate, and /X2, 1^3, ^i, ^2, ^3 and 1)3 are 
constants to be determined. We can take a to be real without loss of generality, uj must 
depend on a if the appearance of secular terms at the third order is to be avoided. Again 
expanding in powers of a, we have 

CU = UJo + Q^UJ2 + . . . , (11) 

where ti^o = —^EqIz^ is the linear dispersion relation for /iq — ^ 00 and uj^ = uj2{k). 

Our expansion for (p satisfies the Laplace equation. Since we are seeking a periodic wave 
train, the boundary conditions @ and (H) do not apply. Inserting our expansions (H), (p^Of ) 
and ([TI|) into Eqs. and (|^), setting a = 1 — 2c, and working to third order in a, we 
obtain 

002 = f3Eo (1 + ^« - k\ (12) 

/i2 = ^ak, (13) 

fi3 = ^a{2a - 1)P , (14) 

ui = -iEo, (15) 

P2 = i{l- ^a) Eok, (16) 

and 



2 V 2 
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The nonlinear dispersion relation is 
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z>3 = h{3a - l)Eok\ (18) 



u = -pEok' 



(19) 



Vp = ^ = -f3Eok i^l - -k'a' + ...] (20) 



For the isotropic case a = 1, the magnitude of the phase velocity 

^ ^„ , 3,22 



is a decreasing function of a when a is small. In contrast, the phase velocity of small 
amplitude gravity waves on deep water increases with amplitude. For a given amplitude 2a, 
the phase velocity is smallest for a = 3/4, or, equivalently, for c = 1/8. 
The form of the periodic wave train is given by 

( = 2a cos{kx — ut) + aka^ cos 2{kx — uot) + -Q;(2a; — l)fc^a^ cos ?){kx — out) + . . . (21) 

To third order in a, the wave is symmetric about the crests, even though the applied electric 

field breaks the right-left symmetry. 

IV. Slow Modulation of Finite Amplitude Periodic Wave Trains 

We will now relax the simplifying assumptions of Section III. Suppose ho is finite and 

that the amplitude is a slowly varying function of position. We set 

oo 

C= E Cni?" (22) 
n=— oo 

and 

oo 

</>= E <^„i?^ (23) 

n=— oo 

where E = e'^i^^-^ot) g^j^^ -g gjyg^ ^j^g linear dispersion relation for films of finite 
thickness, Eq. (^. Since C and are real, C-n = C 'P-n = 0n 

Because the amplitude of the wave is of order e, the coefficients Cn and 0^ are of order 
e". We may therefore write 

oo 

Cn = E e'Cn, (24) 
j=n 

and 

oo 

0n = Ee'>-r (25) 

j=n 

Here n > and the Cn/s and 0nj's do not depend on e. Note as well that Coo = 0oo = and 
that the (^oj's and 0oj's are all real. 

To find a steady, periodic wave train (the analog of the Stokes wave), we would take the 
amplitudes Cnj to be constants and the 0nj's to depend on z alone. We will go further, and 
allow the Cnj's to vary slowly in space and time as viewed from the frame of reference moving 
with the group velocity. As a result, we will be able to find the steady, periodic wave train 
and investigate its stability as well. (Our analysis closely parallels Davey and Stewartson's 
treatment of modulated gravity waves p4|-) To be explicit, we set 



Cn,=Cn,(^,r) (26) 



and 



4>nj = 4>nj{^,Z,T), 



(27) 



where ^ = e(x — Vgt) and r = e^t. The method of muhiple scales will be employed, i.e., x, 



^, t, and T will be treated as independent variables. Using the chain rule, we see that 



d_ 

dx 



d_ d_ 

dx 



and 



d 
It 



d 



(28) 



(29) 



We now insert our expansions for d/dx and into the Laplace equation (p. This results 
in a series of ordinary differential equations for the 0nj's. Applying the boundary condition 
fH), we find that 



■>22 



-- A 
F 



cosh kz' 
cosh kho ' 
cosh 2kz' 



(30) 
(31) 
(32) 

(33) 

cosh/c/iQ coshfc/io coshfc/io' 

where A, F, G and H depend only on ^ and r. It also follows that (pQi and 002 are functions 

of ^ and r alone, while 



and 



-'13 



cosh 2kho 
^ cosh fcz' sinh fcz' 

cosh kho ^ cosh /c/iq ' 



1 , , cosh kz' ^ ,sinh/c2;' ^^coshfc^;' 
--A^i[z) — — iGf:z' — + a- 



^03 



,•9^01 



(34) 



dz d^^ 

The next step is quite laborious but is nonetheless straightforward, and so I will omit the 
details. Inserting our expansions for d/dx, d/dt, ( and into Eqs. and and working 
to second order in e, we find that 

Coi = 0, (35) 
cr 



C: 
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^0 



,12 



• ^ , ho 
i-^G+—Ai., 



,22 



a a 



A' 



(36) 
(37) 
(38) 
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and 



ik 



Equating the coefficient of e^i?° to zero in Eqs. (@) and yields two linear equations 
for ^02 and d(f)oi/dC,. Solving these for ^02, we obtain 



02 



ak\ aahnk — 2 , 
' \A\. 



(40) 



When we equate the coefficient of e^E^ to zero in Eqs. @) and (P), we obtain two 
equations that are consistent only if 



iAr = TA^^ + A\A\^A, 



(41) 



where the constants F and A are given by 

1 d^ujo f3Ef 



2 dk^ 



a 



and 



A 



pk^ 



2 /5 2 \ 2 

— — 1 — a + -a — a cr 
(t2 V2 



1 + 2ahok 



n cr \ 2valvn — aahnk 

1 - 0-^ + \ 9/ U 



(42) 



(43) 



Eq. (^Tj) is the nonlinear Schrodinger equation, and describes the time development of 
the amplitude modulation. It has nonlinear plane wave solutions 



(44) 



where Aq and /t are constants and Vl = F/t^ — A|y4op. Let us recast the solution with k = in 
terms of the original, physical variables. For k = 0, Eq. (|37| ) has the solution C12 = C = 0, 
and we will specialize to this case. We set aeAo = —iEoa, where a is a real constant of order 
e. To second order in a, 



f k\ aah^k — 2 ^ 2A; 
C = 2a cos(kx — ut) + \ — \ ; a H - 



1 - I 1 - -a I 0- 



where 



' E \ 

LU = luq + A { — j a^. 



cos2{kx — cut), (45) 



(46) 



Eq. (^51) describes a periodic wave train of finite amplitude propagating over a metal or 
semiconducting layer of finite thickness. Eq. (|46|) is the nonlinear dispersion relation. 
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In the limit h^k ^ 1, Eq. pq) becomes 



( = 2a cos{kx — ujt) + ako? cos2{kx — ut) + -aka^ + . . . 



(47) 



and Eq. (|46|) reduces to Eq. (|1^). Eq. (^Tf ) reduces to Eq. (|21]) after a trivial vertical shift 
of the origin, and so the results of this section have the correct thick film limit. 
For /lofc < 1, Eq. (||) is 

2 



where 



C = 2a cosikx — cut) — „ [1 — cos 2(/cx — ut)] + . . . 

h^k^ 



u; = vok (l + ^hlk' + j^a' + ... 



(48) 



(49) 



Eqs. (^) and (^Of ) are valid for fca <C {khoY ^ 1- It has been shown elsewhere that the 
equation of motion for small amphtude, long waves is the Korteweg-de Vries (KdV) equation 

1 



3 ho 



(50) 



The KdV equation has finite amplitude, periodic solutions (cnoidal waves) of the form 



( = -2a + 4a cn^ [vr ^K{m){kx - ujt)\m\, 



(51) 



where cn is a Jacobian elliptic function, K{m) is the complete elliptic integral of the first 
kind. 



4a / 2 
Sho Km 



LU = Vok 

and m is given implicitly by the following relation: 



2TT^a 



hlk^' 



(52) 



(53) 



Eqs. (^) - (^) give a good approximate solution to the original equations of motion in the 
limit that a/h^ and h^k tend to zero with a/h^k^ remaining finite as the limit is taken. If 
the amphtude of the cnoidal wave is small enough that ka <^ (kho)^ <^ 1, Eqs. (0) and 



(^) reduce to Eqs. (|48|) and (|^). Thus, in the thin film limit, the nonlinear wave ( ^5|) is a 
cnoidal wave of small amplitude. 

The periodic wave train (^i]) is marginally stable if FA < but is unstable if FA > . 
This means that for a given value of a, the periodic wave train is unstable if kho exceeds 
a critical value we shall denote by Xc = Xc(q^)) and is marginally stable for kho < Xc- The 
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a 



FIG. 2: A plot of 1/xc as a function of a. Xc ^ is zero for a > 2 and a < —1/2. 

critical value Xc = 1.609927... for the isotropic case [a = 1) and is infinite for a > 2 
and a < —1/2. Fig. ^ is a plot of l/xc(ct) as a function of a. It is intereresting to note 
that the smallest value of Xc is obtained for a value of a differing just slightly from unity 
(a = 0.972...). 

For kho <^ 1, the product FA is negative, regardless of the value of a. The periodic wave 
train is therefore marginally stable in this limit. This is in accord with the fact that cnoidal 



waves are marginally stable [|26 . 
IV. Conclusions 

In this paper, the electromigration-induced dynamics of small-amplitude disturbances 
on a low index crystal surface were studied. Periodic wave trains of finite amplitude were 
found, as well as their dispersion relation. These wave trains are unstable for kho > Xc{(^), 
and are otherwise marginally stable. The value of the parameter a depends on strength 
of the material anisotropy, and is equal to 1 if the anisotropy is negligible. For a = 1, 
Xc = 1.609927. . .. For a < —1/2 and a > 2, on the other hand, Xc is infinite and the wave 
train is marginally stable no matter what its wavelength. 

The equation of motion for slow modulations of the finite amplitude, periodic wave train 
was shown to be the nonlinear Schrodinger equation. As is well known, this equation is 



exactly solvable for initial conditions that vanish with sufficient speed as |x| oo ||2^. This 
exact solution has an number of consequences for electromigration-induced propagation 
of surface waves. If kh^ > Xc{c(), a localized initial wave packet of arbitrary shape will 
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eventually disintegrate into a number of envelope solitons and a small oscillatory tail. These 
solitons survive collisions with each other with no permanent change except a shift in phase 
and position. 

What are the prospects of an experimental test of the theory developed in this paper? 
Quite recently, the SEM-induced dynamics of metal surfaces of large area have received some 
attention, but unfortunately these experiments were carried out on polycrystalline films, and 



so they do not permit a test of the theory p8|, ^ 



The experimental situation is more promising in the case of semiconductors. A number 
of beautiful experimental studies of the SEM-induced dynamics of single crystal silicon films 
have been carried out, but on initially planar, wcma/ surfaces 0. These studies revealed 
that a vicinal surface is unstable against step-bunching for one current direction, and that 
it is stable for the opposite current direction. Experimental studies of the dynamics of 
perturbed, low-index silicon surfaces that are subject to high electric fields have not yet 
been carried out. 

To test the theory, a singly-periodic periodic wave structure could be etched into a low- 
index silicon surface. After perturbing the surface in this way, a high electrical current 
would be applied parallel to the ripple wavevector. The subsequent surface dynamics could 
be imaged by scanning tunneling microscopy, as has already been done for vicinal surfaces. 

The analog of the Stokes wave would be obtained by etching many parallel ripples into 
the wafer surface before the application of current. In this way, the amplitude dependence 
of the phase velocity could be determined and compared with the theoretical prediction. To 
observe the development of an envelope soliton (or solitons), on the other hand, a few tens 
or hundreds of ripples would be etched into the film initially. This would produce a "wave 
packet" that would ultimately disintegrate into one or more envelope solitons and a small 
oscillatory tail under the action of an applied current. 
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